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Abstract 

Economic indicators time series are usually complex with high frequency data. The tra- 
ditional time series methodology requires at least a preliminary transformation of the data 
to get stationarity. On the other hand, the Robust Bayesian Dynamic Model (RBDM) docs 
not assume a regular pattern and stability of the underlying system but can include points of 
statement breaks. In this paper, we estimate the Consumer Price Index and the Economic 
Activity Index of Puerto Rico using a RBDM with observational and states variances to 
model the outliers and structural breaks in the time series. The results show the model de- 
tects structural changes in both series, some of them just before the beginning of a recession 
period in Puerto Rico's economy. 
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1 Introduction 

In this paper we present a Robust Bayesian Dynamic Model (RBDM) to model the outliers and 
structural breaks of a time series. This model permits to take into account not only the outliers 
and structural breaks of the historical series, but also allow us to do posterior inference using 
a Bayesian approach. A linear trend model, also called linear growth model, with robust priors 
for the distributions of the state and observational precisions errors was fitted to the historical 
series of two economic indexes widely used to describe the economic situation of a country: the 
Consumer Price Index (CPI) and the Economic Activity Index (EAI). The former measures 
inflation through the fluctuations of prices, and the latter measures the real economic activity. 

In Puerto Rico, the CPI is produced and published by the Department of Labor and Human 
Resources (DTRH) and the EAI is produced and published by the Government Development 
Bank (GDB). The CPI represents the average monthly change in prices of the items and services 
that urban families usually buy, and it is used to study the behavior of inflation. The CPI began 
as an index of cost of living in Puerto Rico, published as the Index of Cost of Life for the Families 
of Workers (see Department of labor and human resources (2008)), and was developed in the 
early years of the 1940's as an index of cost of living for working families. This index was reviewed 
approximately every ten years. For the first time, in the Study of Income and Expenses of 1977, 
information was collected about the urban families in Puerto Rico, not only of working families, 
including the self-employed and the pensioners, among others. The use of the original index 
was discontinued in 1980. With the new information collected in 1977, a new index was born, 
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the index of prices to urban consumers, widely known as the Consumer Price Index (CPI). 
After the Study of 1977, a Study of Income and Expenses is done again in the years 1999-2003. 
Since the elapsed time between studies was more than twenty years, in 1990 a few adjustments 
to the items and services included in the basket that makes up the CPI were made, including 
the addition of new products. The continuing revision of the index is extremely important so 
that it reflects accurately the changes in tastes and preferences of the population, as well as 
the inclusion of new products in the market that are not represented in the previous basket. A 
major change in the CPI calculation methodology is implemented in March 2010, incorporating 
the use of the geometric mean for the aggregation of prices, a measure used by the majority 
of countries, including the United States. Currently, the base year for the basket of goods for 
the CPI is the year 2006 (the value of the CPI for 2006 is 100). This basket is composed of 
the following major groups: food and beverages, housing, apparel, transportation, medical care, 
entertainment, education and communication, and other items and services. These groups are 
similar to the ones in the United States basket (see Department of labor and human resources 
(2008)). 

The Economic Activity Index (EAI) of the Government Development Bank, on the other 
hand, is a monthly index that includes the behavior of four economic indicators: total of salaried 
employment (thousands), sale of cement (million bags), fuel consumption (millions of gallons) 
and electricity generation (million KWH). Until March 2012, one of the indicators was electrical 
energy consumption, but it was replaced since then by electricity generation. The Government 
Development Bank for Puerto Rico (2012) attributed this change to "the electricity generation 
data is obtained in a timelier manner, is more reliable because they have fewer revisions, and 
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there is a slightly higher correlation with the Gross National Product (GNP) of Puerto Rico" . 
The Government Development Bank for Puerto Rico (2012) reported that the correlation of this 
new index with the GNP is 0.97, which means that they are highly correlated. 

We consider the monthly series of the CPI from January 1980 to December 2012; the monthly 
data for the AEI goes from January 1980 to December 2011. The CPI and the AEI have had 
different abrupt changes in the last years. Therefore, it is interesting to have a model that allows 
the detection of these changes in both indexes. As the CPI and AEI are two important measures 
of the economy, we will also investigate the relationship between the changes the two series have 
had in Puerto Rico. 

The paper is organized as follows. Section 2 shows the RBMD prior variances specification 
for the RBDM. Section 3 presents the RBDM to modeling the CPI and the AEI. Finally we 
have the conclusions and remarks in Section 4. 

2 Model Specification and modelling outliers and structural breaks 

A univariate Dynamic Linear Model (DLM) is specified (see West k, Harrison (1997)) by the set 
of equations: 

y t = F t e t + u t v t ~N(0,V t ), (1) 
e t = G t 0t- 1 +UH u t ~N(0,W t ), 

where t = 1 : T. The specification is given by the prior distribution for the initial state 0q. 
This is assumed to be normally distributed with mean tjiq and variance Cq. yt and 9t are m 
and n-dimensional random vectors and F t , G t , Vt and W% are real matrices of the appropriate 
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dimension, yt is the value of an univariate time series at time t, while 6t is an unobservable state 
vector. On the other hand, the scaled Beta2 prior for the precision A = 1/r 2 is the following: 



fW = rirl^(i + ^' A>0 (2) 

where small values of (3 are considered in order to have heavy tails priors for robust inference. 
This paper consider the Student-t density coupled with the scaled Beta2 (see Appendix B) for 
modelling the observational and states variances in RBDM (see Fuquene, Perez & Pericchi 
(2013)). Therefore let 6 ~ Student-t(/x, r, v) where v are the degrees of freedom, fi the location 
and r the scale of the Student-t density: 



2^ 



>+l)/2 



,r(^) = -^l + -(_^ j , v > 0,-oo < n < oo,-oo < 9 < oo, (3) 

where fci = ^ t^/^ ■ We have that tt(6>) = J °° vr(^|r 2 )7r(r 2 )(iT 2 therefore 



'/3^/(0- /[ /)9+ 1 /2 2 Fl(p + ^,g + l/2,(t; + l)/2+p + g,l-/3^/(e-/ i ) 2 ) if 6> 7^ >u, 



n(9) 



k fciBe(g + l/2,p + u/2)/Be(p,g) if = p. 

with 2Fl(a, 6, c, z) the hypergeometric function (see 15.1.1 of Abramowitz & Stegun (1970)) 
and we have that ir(6) is the Student-t-Beta(w,p,q,/3) (see Fuquene et al. (2013) for the proof of 
this result). For example ifv=p = q = lwe have the Student- t-Beta2( 1,1,1,/?)) such as: 
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< d) = — t\ — i\* (4) 

We can see in Figures 1 and 2 the student-t-Beta(l,l,l,/3) prior has heavier tails than the Cauchy 
prior. That means these priors can detect outliers even far from to these obtained with student 
priors. Also, using these priors a simple Gibbs sampler can be used because all full conditional 
densities in the gamma hierarchical parameters have gamma distributions. 



to 




e 



Figure 1: Comparison of the Student-t-Beta2(l,l,l,l) ,Cauchy(0,l), Normal(0,2.19) priors. 
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Figure 2: Comparison of the tails of the Student-t-Beta2(l, 1,1,1) ,Cauchy(0,l), Normal(0,2.19) 
priors. 



In order to model the CPI we use the Student-t-Beta(u,q,p,-g) (using the Beta2 prior for the 
precision A = 1/t 2 ). Wt t i denotes the ith diagonal element of Wa, i = 1, ■■■,n the hierarchical 
Student-t-Beta(u,q,p,i) prior can be summarized such as: 



V t 1 =XyUJy : t, 

X y \q ~ Gamma(g, (Ppy)^ 1 ), 
^y,t ~ Gamma(f /2, 2/v), 
p y ~ Gamma(j>, 1), 



M,i\q ~ Gamma(c7, (/3p eA ) _1 ), 
LjJ e,t l ~ Gamma(u/2, 2/v), 
pe,u ~ Gamma(p, 1), 



For each t, the posterior distribution of uj Vy t (i.e. <jJe,u) contains the information of outliers 
and abrupt changes in the states. Values of u Vj t (i.e. ^e,u) smaller than one indicate possible 
outliers or abrupt changes in the states (See Petris, Petrone & Campagnoli (2010)). A Gibbs 
sampler is implemented using the posterior distribution of parameter and states of the model 
specified above. 

We have the annual Consumer Price Index in Puerto Rico on a log-scale in Figure 3. 




i i i i i 

1985 1990 1995 2000 2005 2010 
Time 

Figure 3: Annual Consumer Price Index in Puerto Rico on a log-scale, 1984-2010. 

The natural choice for modelling the logarithm of the CPI using a DLM is a local linear trend 
model or also called linear growth model which fit the trend and slope of the CPI logarithm. 
The linear growth model is the following: 
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y t = fit + v t , v> t ~N(0,V t ), 

fit = M-i + 6-1 + w tj i, U) t ,i ~ A^(0,cr t 2 j(U ), (5) 

6 = 6-1 + ^t,2, w t)2 ~ iV(0, of^), 

with uncorrelated errors i/j, u^i and 0^2- This is a DLM with: 
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3 The Robust Bayesian Dynamic model for the CPI 



In this section we implement the proposed model. First a motivating example is showed where 
the proposed model is illustrate using the annual logarithm of the CPI. Next we apply our 
approach to the CPI and AEI monthly indexes. 

3.1 A Robust Bayesian Dynamic model for the annual logarithm of the CPI 

We implement the proposed model for modelling the annual logarithm of the CPI. We use a 
Student-t with four degrees of freedom and a non-informative Gamma for modelling the outliers 
and changes in the variance states. The choice v = 4 degrees of freedom is not new. For example 
Gelman (2004) recommend the Student-t prior with four degrees of freedom in order to obtain 
robustness in Bayesian applications. 
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Figure 4: Outliers and structural breaks in the annual logarithm Consumer Price Index in Puerto 
Rico using the robust approach 
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In the proposed model we use a Student-t-Beta2 where p = q = 1, and 1//3 = 10000 for the 
posterior robustness inference. We obtained convergence of all parameters using 30000 iterations 
after a burn- in period of 10000 iterations. The Figure 4 displays the proposed model we can 
see there is an outlier at 2000. The trend shows changes in 1990 and 2001 and the slope has a 
change in 2003. The proposed model has taken into account not only the outliers but also the 
structural breaks for the annual logarithm of the CPI. 

3.2 Bayesian Dynamic Robust Model for the CPI and EAI monthly indexes 

The motivating example of the previous section shows how the proposed model works. However, 
the most interesting case is when the model is fitted to the monthly series of the two indexes 
studied. Since the CPI and AEI indexes could be related, we apply the proposed model to both 
series, in order to have a broader analysis of the changes during the period studied. 

The EAI is one of the indicators used to detect a recession period in the economy, since it is 
highly correlated to the Gross National Product (GDP). When the GDP decreases for more than 
six months consecutively, the economy could be going through a recession period. In Puerto 
Rico, the recession periods identified during the last 35 years are: 1980-82, 1990-91, 2001-02 and 
since 2006. 

Figure 5 displays the results using the proposed model for the CPI for the period of January 
1980 to December 2012. Important remarks can be made. The residuals in the bottom of Figure 
6 are given by h = Vt — E{F9t\yi-.T)- By looking at the residuals it can be seen there is only 
one outlier in January 2009. The slope is dynamically changing with only a sudden jump in 
September 1981, within the first recession period. 

The trend has different jumps, the most dramatic one in September 2005 with E(uo,t 1 \yi-.T) = 
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Table 1: Posterior mean of we,t x for the monthly logarithm of the CPI. 



month/year E(uig^ \Ui-.t) - CPI 



May/1980 


0.26512911 


July/1989 


0.15719286 


September/ 1990 


0.25577646 


April/2005 


0.20028956 


September/2005 


0.07095877 


December/2005 


0.26580040 


April/2006 


0.25441340 


October/2006 


0.13192695 


December/2006 


0.12591796 


June/2008 


0.26553021 


November/2008 


0.11320069 


December/2008 


0.13743595 



0.07 and some other abrupt changes in the precedent years. Even though the Consumer Price 
Index is an indicator of inflation and not recession, this dramatic change could have been an 
"alarm" for the economic recession Puerto Rico is facing since 2006. Other dramatic changes 
are found in May 1980, July 1989 and September 1990, all of them just before the beginning of 
a recession period in Puerto Rico, as shown in the previous paragraph. The structural break of 
1980 could be also tied to the changes in the methodology to compute the CPI, as discussed in 
the Introduction. The change in the trend of the CPI at the years 1989 and 1990 could also be 
related to the implementation of the "Joint Committee on Taxation" in the United States, and 
its effect on the island. 
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Table 2: Posterior mean of uie^i for the monthly logarithm of the EAI. 



month/year 


E(u g , tl \vi-.t) - EAI 


July/1980 


0.4879160 


March/1983 


0.3089463 


December/1987 


0.4912573 


September/1989 


0.2574237 


December/1989 


0.2864106 


January/1990 


0.4640351 


July/1996 


0.2484716 


Decembcr/1996 


0.3511382 


September/1998 


0.1563733 


October/1998 


0.3810037 


November/1998 


0.2929802 


December/1998 


0.2276613 


December/2001 


0.4530024 


September/2005 


0.3313795 



Figure 6 shows that the abrupt changes in the trend for the EAI are much related with those 
found in the CPI. In particular, the most dramatic ones are presented in September 1989 and 
December 1989, just before the beginning of the recession period of 1990, July 1996, September 
1998, December 1998, and September 2005, coinciding with the break in the CPI series, just 
before the beginning of the 2006 recession. Figure 7 displays the relationship between the indexes 
for the periods 2000-2005 and 2005-2010. In particular, the changes in the period 2000-2005 
for the two time series are very similar. The abrupt changes in both indexes during the second 
period may be the consequence of the economic crisis that Puerto Rico has been suffering since 
2006. 
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In conclusion, the RBD models fitted to the series of the Economic Activity Index and the 
Consumer Price Index of Puerto Rico, during the period 1980 - 2012, detect abrupt changes in 
both series that may be tied to the beginning of a recession period in the economy. 
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Figure 5: Outliers and structural breaks in the logarithm monthly Consumer Price Index in 
Puerto Rico 
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Figure 6: Outliers and structural breaks in the logarithm monthly Economic Activity Index in 
Puerto Rico using the usual approach 
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Figure 7: Comparison structural breaks for the monthly Consumer Price and Economic Indexes 
in Puerto Rico 
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4 Conclusions 

In this paper we present a robust Bayesian dynamic model for the CPI and the EAI. We fitted 
a linear trend model with robust priors for the distributions of the state and observational 
precisions errors. The Bayesian dynamic robust model has the quality to detect outliers and 
historical structural breaks. The changes of the CPI trends are associated with those changes 
obtained for the EAI series. In fact, the structural changes in both series have a contextual 
historical and economical meaning. Some of these structural changes in both series coincide 
with the beginning of a recession period in the economy of Puerto Rico. Finally, the proposed 
model has the feature that it produces credible intervals that are not constant over time. 
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A Details for the Markov Chain Monte Carlo algorithm 

This appendix presents the MCMC scheme we use in the paper. 

• FFBS algorithm for the states. 

In order to obtain posterior inference on the state parameters 6t in the model (5), we use 
the forward filtering backward sampling (FFBS) given in Fruwirth-Schnatter (1994) which 
is practically a simulation of the smoothing recursions. The FFBS works for the model 
(5) as follow: 

1. Use the Kalman Filter equations for (5). Let mo and Co (known) with (0o\Do) ~ 
N(mo, Co) and 



9t\yi-.t-i ~ N(m t -i, C t -i) 



(6) 



— The one step predictive distribution of 9t given yiu-i is Gaussian (6 t \D t -i) 



N(at,Rt) with parameters: 



a t = G t m t -i; 




(7) 
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— The one step predictive distribution of yt given yi-.t-i is Gaussian (yt\Dt-i) ~ 
N(ft,Qt) with parameters: 

ft = F t a t ; Qt = F t R t F t + V t . 

- The filtering distribution of 9t given yv.t-i is Gaussian (9t\Dt) ~ N(m,t, Ct) with 
parameters: 

m t = at + A t e t ; C t = R t - A t Q t X t (8) 

where A t = R t F t Q^ 1 , and e t = y t - f t . 

2. At time i = T sample 6> T from iV(0 T |mt, C t ). 

3. For t = (T - 1) : sample 9 t from JV((9 t |mJ, C t *)with 

m* t =m t + B t (9 t+l - a t+1 ) C* t =C t - B t R t+1 B' t 

where B t = C t G' t+1 R^ v 

Full conditionals for the observational and state variances in 
the DLM. 



T 1 \ / T 1 

A y |... ~Gamma ( q + -, -SSy* + f3p y J , A^|... ~Gamma U + -, -SS* ei + Pp e ,t, 



Wy t t\... ~Gamma I — - — , 1 u e>ti \... ~Gamma I — - — , - 

p y \... ~Gamma(p + q,f3X y + 1) , P0,u\— ~Gamma (p + q, /3Xg ti + 1) 



where SSy* = ELi u vAvt ~ F t^t) 2 and SS* e - = £f =1 w* )ti (0* - (G^),) 2 for i = 1,2 
and t = 1....T. 
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B Scaled Beta 2 density 

A Beta2 random variable is equal in distribution to the ratio of two gamma-distributed 
random variables having shape parameters p and q and common scale parameter (3: 



t 2 ~ Gamma(p,/3/p) (9) 
p ~ Gamma(g, 1) (10) 



where Gamma(a, b) denotes the Gamma distribution: 



p(x\a,b) = —\- a x a - 1 exp{-x/P} a>0,b>0, (11) 
1 (ajp 

with f3 the scale parameter. The scaled Beta2 prior can be defined such as: 



... r(p + g) i V/?y _ v n 

f(T) = W/ 1 + ^'' T>0 - (12) 
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